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We analyze the prospects for observing the light and heavy CP-even neutral Higgs bosons (</> = 
h and if) in their decays into flavor violating bs channels (including charge conjugation) at the 
proposed Large Hadron electron Collider (LHeC), with ^/s « 1.3 TeV, in the framework of a 2- 
Higgs Doublet Model (2HDM) Type-III, assuming a four-zero texture in the Yukawa matrices and 
a general Higgs potential. We consider theoretically consistent scenarios in agreement with current 
experimental data from flavor and Higgs physics. We investigate the charged current production 
process v e <fq in presence of flavor violating decays of the Higgs bosons, that lead to a 3-jets + fir 
signature. We demand exactly two jets, one tagged 5-jet and one light-flavor jet, all in the central 
rapidity region. The remaining jet (originated by the remnant quark q) is tagged in the forward 
or backward regions and this together with a central jet veto (not more than one light-flavor jet) 
are essential criterions to enhance the signal-to-background rates. We consider the most relevant 
Standard Model (SM) backgrounds, treating c-jets separately from light-flavor and gluon ones, while 
allowing for mis-tagging. We find that the SM-like Higgs boson, h , would be accessible within several 
parameter configurations of our model at approximately the l-2cr level with 100 fb -1 of data. We 
also find that the heaviest neutral Higgs boson, if, with mass up to 150 GeV, would have a ler 
significance for the same data sample. At the end of the LHeC running, one would have ten times 
data accumulation and for all the Higgs masses the significances are enhanced so as to allow for 
detection of both the h and if state. Hence, one of the most viable extensions of 2HDMs with Flavor 
Changing Neutral Currents (FCNCs) generated at tree level but controlled by a four-zero texture 
approach in the Yukawa matrices, as opposed to the adoption of ad hoc discrete symmetries, could 
be put under scrutiny at a future ep machine. 
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I. INTRODUCTION 


The Standard Model (SM) is well established by now after the discovery of a Higgs boson by the ATLAS [1] and 
CMS 0 experiments at the Large Hadron Collider (LHC). However, when one considers some theoretical aspects of 
the SM, for example, lepton number violation, which is already manifested in the form of small but non-zero neutrino 
masses and actively searched for in the other two fermionic sectors (e.g., in (i -A ey decays by the MEG experiment 
Q and in H-physics by BaBar [4] and Belle IB one necessarily has to postulate New Physics (NP) Beyond the SM 
(BSM). By combining the evidence of Higgs states and the presence in Nature of flavor violation, then one is well 
motivated in searching for evidence of BSM physics in a context where the two aspects merge, i.e., in flavor violating 
Higgs boson decays. In general, limited to the Higgs sector, several BSM scenarios have been invoked by introducing 
extra Higgs singlets, doublets and/or triplets. As the Higgs boson discovery at the LHC is consistent with a doublet 
structure, we refrain here from considering BSM constructions with either of the other two aforementioned Higgs 
representations. Therefore, in staying with multi-Higgs doublet structures, the simplest of such Higgs scenarios is the 
so called 2-Higgs Doublet Model ( 2 HDM) [ 6 , 7i], which will be the theoretical focus of our study. 

Among the many phenomenological sides of a 2 HDM, we are indeed concerned here with flavor violating Higgs 
boson decays in the quark sector, building on previous works of ours, see, e.g., Refs. Bum. However, plentiful of 
studies, also including lepton flavor violating scenarios, exist, some specific to 2HDMs and others adaptable to their 
case: see for an incomplete list, e.g., Refs. |16l-l32]. 

The actual search scope for Higgs bosons in flavor violating modes at the LHC has also been actively studied, see 
It has also been investigated in the context of a future e + e - [I^ and 77 collider [l2|. Prospects at a future 
hadron collider have been investigated in [38| . 

Herein, we are particularly motivated by a possible enhancement of flavor violating quark decays (0 -A bs) of 
intermediate mass Higgs bosons (below the top-quark mass) and we will focus on the possibility to access such 
signatures at the possibly upcoming Large Hadron electron Collider (LHeC). The LHeC facility (44| presently discussed 
as possible at CERN in the near future would be a Deep Inelastic Scattering (DIS) experiment at the TeV scale, with 
center-of-mass energy of around 1.3 TeV. In comparison to the another recently closed (in 2007) DIS experiment (the 
Hadron-Electron Ring Accelerator (HERA) [45| at DESY of around 320 GeV in energy with integrated luminosity 
of around 0.5 fb _1 ), the LHeC might deliver data samples of approximately 100 fb -1 and with a higher detector 
coverage. Further, the overall kinematic range (in x and Q 1 2 ) accessible at LHeC is 20 times larger than at HERA. 

While the primary task of a collider like the LHeC will be in-depth studies of QCD, the machine also affords some 
scope to study Higgs bosons decaying via flavor violating processes [46). Our objective in this paper is to study the 
feasibility of finding two CP-even neutral Higgs bosons, the SM-like Higgs state h and its heavier counterpart iL, at 
the upcoming LHeC assuming as BSM framework a 2 HDM Type-III (henceforth 2 HDM-III for short), which embodies 
a four-zero texture approach in the Yukawa matrices as the mechanism to control Flavor Changing Neutral Currents 
(FCNCs). As we have shown in previous analyses, this is precisely the framework which establishes 0 -A bs + c.c. as 
an hallmark signature of quark flavor violation in the Higgs sector whose detectability is under investigation at the 
LHC and future e + e - and 77 collider. 

The plan of our paper is as follows. In the next Section we describe briefly the theoretical structure of the 
2HDM-III with a four-zero texture embedded in the Yukawa matrices. In Section III, we demarcate the allowed 
2HDM-III parameter space in presence of both theoretical and experimental constraints. In Section IV, we explain 
the characteristics of the Higgs boson signal from charged current production. We introduce the most important SM 
backgrounds and finally, we carry out both a parameter space scan and a signal-to-backgrounds analysis for some 
characteristic benchmarks by adopting a simple cut-based optimization to isolate 0 -A bs + c.c. events. In section V, 
we recap and present our conclusions. 


II. THE HIGGS-YUKAWA SECTOR OF THE 2HDM TYPE-III 

In the 2HDM, the two Higgs scalar doublets, = (0j~,05*) and = ( 02 ,02*)? have the same hypercharge +1 
such that both couple to the same quark flavor. Since a specific four-zero texture is implemented as a flavor symmetry 
in the Yukawa sector, this is the mechanism which controls FCNCs so that discrete symmetries in the Higgs potential 
are not needed |8l4l2). Then the most general SU(2)lxU (l)y invariant scalar potential, following Q, can be written 


1 In fact, also top-quark flavor violating decays into charm quarks and Higgs bosons are currently under investigation at the LHC. 

2 Current experimental results at the LHC include both ATLAS mi and CMS /l2] analyses. 
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as: 

V{*i, $ 2 ) = + H.C.) + (1) 

+iA 2 ($!$ 2 ) 2 + A 3 ($I$i)(^ $ 2 ) + A 4 ($ t 1 $ 2 )(^$i) 

+ Qa 5 ($1$ 2 ) 2 + (a 6 ($1$!) + A 7 ($t$ 2 )) ($1$ 2 ) + H.c.^ , 

where all the parameters are assumed to be rea@, including the Vacuum Expectation Values (VEVs) of the scalar 
fields, hence there is no CP-violation. In general, by introducing some discrete symmetry <Li —>■ and <£>2 —> —$2 , 
the scalar potential does not have the contributions of A 6 and A 7 . 

It has long been proposed that there are four possibilities to satisfy the Paschos-Glashow-Weinberg theorem m in 
2HDMs [ 6 ],[7]. These are defined as follows: Type I (where one Higgs doublet couples to all fermions); Type-II (where 
one Higgs doublet couples to up-type quarks and the other to down-type ones); Type-X (also called “Lepton-specific”, 
where the quark couplings are Type-I and the lepton ones Type-II); Type Y (also called “Flipped” model, where 
the quark couplings are Type II and the lepton ones Type-I). With these two scalar doublets, there are eight fields 
but only five of them are physical (pseudo)scalar (“Higgs”) fields, which correspond to: two neutral CP-even bosons 
h (lighter) /H (heavier), one neutral CP-odd boson A and two charged bosons H ± . The mixing angle a of the two 
neutral CP -even bosons h and H is another parameter of the 2 HDM model. In total, the 2 HDM model can be 
described by the parameters <a, /3 (where tan [3 is the ratio of the VEVs of the two Higgs doublet fields) and the 
masses of the five Higgs particles. With these inputs one can estimate all the parameters that are present in the scalar 
potential, to be specific, the A’s. These A’s (together with various scalar mass parameters) enter the expressions of 
the theoretical constraints like: vacuum stability, unitarity, perturbativity and also various EW Precision Observables 
(EWPOs), for example, the oblique parameters. All these 2 HDM types are fully compatible with the SM-like Higgs 
boson discovery. 

The flavor sector of 2HDMs is testable in low as well as in high energy collider experiments. The tests have been 
carried out in the most general version of a 2HDM with a Yukawa four-zero texture, wherein the Yukawa couplings 
are proportional to the geometric mean of two fermions masses, gij oc yj(mirrij)xij Bill. As it was mentioned, a 
consequence of this is that the terms of the scalar potential including Ac and A 7 should now be taken into account. 
This leads to tri-linear and quartic self-couplings of the scalar fields [ 8 |, [l2j affecting the model phenomenology in 
one loop processes via di-Higgs and tri-Higgs topologies, both in production and decay processes. It has been shown 
that the EWPO p can deviate from experimental bounds at one loop level, as long as the mass difference between 
charged Higgs bosons with CP-even/CP-odd masses is large, irrespective of the value of A 6 and A 7 . Hence, some level 
of degeneracy between one neutral and the charged Higgs states is a precondition on the 2HDM spectra . In our 
construction, the Yukawa Lagrangian [ 8 ] is given by: 

C Y = - (y^Ql^ur + Y 2 u Q l $ 2 u r + VTMidrt + Y 2 d Q L $ 2 d R + yRl^Ir + Y^LM^j , (2) 


where ^ 1,2 = (</> 1 2 ? ( t > i 2 ) T re f er to the two Higgs doublets, ^ 1,2 = ^ 2^1 2 - Besides, the fermion mass matrices after 
EW symmetry breaking are given, from Eq. ©, by: Mf = -^=(viY/ + V 2 V/), / = fo d, assuming that both 

Yukawa matrices Yl and Y<[ have the four-zero-texture form and are Hermitian [10, 0, Hij. After diagonalisation, 
Mf = Vj L MfVfR , one has Mf = ^(^ 1 ^/ + ^ 2 Y/), where Y/ = Vj L Y/VfR. One can obtain a compact and generic 
form for the rotated matrix Y$ 


Y f 


m{rrij 


m{rrij 


i'dl 


( 3 ) 


JZJ V V 1 

where the y’s are unknown dimensionless parameters of the model. Following one has a generic expression for 
the couplings of the Higgs bosons to the fermions given as 




—Ui (rn dj XijP R + m Ui YijP L ) dj H+ + YH^L Zij 17ll R H + + H.c. 




( 4 ) 


3 The n \ 2 5 A5, A6 and A7 parameters could be complex in general, but for simplicity we assume these parameters to be real. 

4 We have shown in several parametrisations that this structure corresponds, as a particular case, to the Cheng and Sher ansatz HOI 

ill. 
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where (/>{• (</> = h , if, A), Yij and Zij are defined as follows: 


h f ^ 5ij+ 1^,-A ui 


H?j = tUy 


V2f(X) 

(d+*dr) 

V2f(X) 


4 = hfj(d -+e,X^Z), 
nil = H?j(d -»• e,x Z), 


= -XSij + 


/(X) , 


a /2 V 


—~Xijt Aij — Afj(d -t£,X—>- Z), A™j — A t j(d —► u, X —► X), 


Ki=a5n - 


(Ch + YQt) J. n „ u _ £U , _ ^ (ft - ^1) 


V^/(V) 


^ ^=^ + ' V 2 / ( y) 


Xij 5 


*y = E( Vc: 

1=1 

** = t 


m dl c /(X) lm dl „ d 


X-hj- r- m 

m dj V2 V rn d , 


Xij 


Y Su — 


f(X) f^i 


■Xu 


Zij = 


m h 

mu 


y/2 V rn Ui 

f(Z) [m~ i 


(Vckm)zj, 


- /— 4 / 

>/2 V mi j 


Xij 


(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

( 11 ) 


where /(x) = y/l + x 2 , ^ are related to the trigonometric ratios (i.e., cos <a/sin/?, sin <a/sin/?, cos <a/cos/?, 
sin a/ cos/?) and the parameters X, Y and Z can be related to tan/? or cot/?, according to the various incarna¬ 
tions of 2HDMs 0 (see the Table HJ). Taking into account that the Higgs-fermion-fermion (0//) coupling in the 
2HDM-III is written as ^hdm-HI = #2 HDM-any + where ^HDM-any is the coupling </>f f in some of the 2HDMs 
with discrete symmetry and Ag is the contribution of the four-zero textural, it was pointed out in [10] that this La- 
grangian could also represent a Multi-Higgs Doublet Model (MHDM) or an Aligned 2HDM (A2HDM) with additional 
flavor physics in the Yukawa matrices. 


2HDM-III 

X 

Y 

Z 

a 

d 

d 

Ch 

dn 

dn 

2HDM-I-Like 

-ctp 

ctp 

-ctp 

C(xj Sp 

Col/ Sp 

Col/ Sp 

Sol/sp 

Sol/sp 

Sol/sp 

2HDM-II-Like 

tp 

ctp 

tp 

Col/ Sp 

Sol/ Cp 

Sol/ Cp 

Sol/sp 

Col/ Cp 

Col/ Cp 

2HDM-X-Like 

-ctp 

ctp 

tp 

Cql/ Sp 

Col/ Sp 

Sol / Cp 

Sol/sp 

Soc/Sp 

Col/ Cp 

2HDM-Y-Like 

tp 

ctp 

-ctp 

Col / Sp 

Sol / Cp 

Col/ Sp 

Sol/sp 

Col/ Cp 

Sol/sp 


TABLE I. Parameters (</> = h, H , A with / = u, d, £ ), A, Y and Z as defined in Eqs. (I5l)- (HT1) for four versions of 
the 2HDM-III. These parameters are related to the interactions of Higgs bosons with the fermions given in Eq. Here, 

tp = tan /?, ctp — cot /?, Sa - sin a, c a = cos a, sp = sin /?, cp = cos /?. 

Here, we consider three different incarnations of the 2HDM-III, which correspond to to the four 2HDM types already 
described except the lepton specific one, as here leptonic Branching Ratios (BRs) are dominant, whereas we intend 
to look for an enhancement in the Higgs to bs decay because of flavor violation. We will finally show that, in different 
scenarios of the 2HDM-III, a substantial enhancement of the BR((/> -4 bs) (including charge conjugation) is possible. 
We do so first via a parameter scan of the 2HDM-III at the inclusive level, followed by the detailed event generation 
analysis of some benchmark scenarios amenable to phenomenological investigation. 


III. THE 2HDM-III PARAMETERS AND BENCHMARKS 

In this section, we will perform a parameter scan of the 2HDM-III of interest from which we will extract our 
benchmark scenarios, all of which will be studied in our final signal-to-background simulations, albeit we will show 
detailed results only for a subset of these for reasons of space. 

First, we ought to explain the constraints we have enforced in our analysis. As for the experimental ones, we have 
taken into account recent experimental bounds from flavor physics [9|, [lOj: i.e., from B -A rz/ r , D —/iz/, D s -A 


5 For example, one can recovers the Yukawa interactions given in Refs. f50l452ll with x{j — 0- 
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the semileptonic transition B —>• Dtv t , the inclusive decay B L?o — Bo mixing, B s /i + (i~ and the radiative 

decay Z —>> bb. (We have also imposed EWPO limits.) On the theoretical side, we have enforced perturbativity, 
triviality, vacuum stability and unitarity constraints ||f53|. In all the constraints mentioned above, the charged Higgs 
bosons masses are the crucial parameter, as diagrams with B[ ± states co-exist alongside those involving the SM W± 
exchange diagrams. In this connection, alongside flavor and EWPO constraints, we have also accounted for those 
stemming from Tevatron and LHC searches [10, EH, Edl loci pl. 


A. Parameter scan 

We scan the parameters space of the model and we only consider as viable the points that avoid the aforementioned 
theoretical and experimental bounds and that are fully consistent with the most recent results of Higgs physics from 
LHC. Taking into account that our model provides interesting new physics in the form of a substantial enhancement 
of the decay </> sb + h.c. (</> = h, iL), as a direct consequence of the off-diagonal terms of the texture of the Yukawa 
matrices. In our scanning, we ask that BR(0 —» sb+ h.c.) ~ 0.01 to 0.1, keeping the decay (j) bb dominant. Firstly, 
we scan the off-diagonal terms of the Yukawa matrices and after we chose some interesting sets of the y’s parameters, 
which are consistent with the flavor physics constraints and Higgs physics bounds used in the analysis of |9l-[lll| , where 
we have shown that several meson-physics processes are very sensitive to charged Higgs boson exchange, and the 
off-diagonal terms of Yukawa matrices given in the Eq. © are kept constrained in the following range: 

- 0.06 < (Xn )23 < 0.3, -0.3 < (Xnhs < 0.5. (12) 

Secondly, we fix the y’s parameters and the masses of the following Higgs bosons, rrih = 125 GeV, uia = 100 GeV 
and m H ± = 110 GeV. We run the mass ra# of the Higgs boson H from 130 GeV up to 200 GeV. Therefore, we can 
reduce the study of the parameter space to that of the couplings X and Y only, which are constrained strongly by 
the inclusive radiative decay B X s j through the following bound: 

- ij £ £ a3 ' <i3) 

where V 33 , Y 32 are defined in Eqs. m-m and Vtb and V ts are elements of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix. From the constraint in Eq. (USD, we can define the allowed region for two general cases, (a) For the case I 
defined by: X = —Y or X = V, 0.1 < cos (/3 — a) < 0.5, and fixing the parameters of Yukawa matrices, Xkk = 1-5 

{k= 2,3), X22 = !- 8 > X33 = 1 - 2 , X23 = 0-2, X22 = °-5, X33 = L2 > X23 = 0.1. One can see in Fig. H] in the left panel the 
allowed region for Y = —X <15 and for the case X = Y <20. This region could represent the case of the 2HDM-I 
plus a deviation given for the flavor symmetry of the Yuwaka matrices, (b) For the case II given by: X » Y, with 
cos(/? - a) = 0.1, X22 = 0-5, X33 = !- 4 > X22 = 2 > X33 = 1 - 3 > X23 = -0-53, X23 = °- 2 > X22 = 0.4, X33 = L2 > we can 
see that the large values for parameter X is permitted. This case could be the incarnation of both the 2HDM-II and 
2 HDM-Y (or flipped model) plus a deviation given by the four-zero texture of the Yukawa matrices. Considering 
these criteria, we chose three interesting scenarios from the versions of 2HDM-III given in the Tab. HI Scenario lb 
which is related to the 2HDM-I-Like, with cos (/ 3 ) =0.5, Scenario Ila is the case 2HDM-II-like, with cos (/ 3 ) =0.1, 
and Scenario Ya is the case 2HDM-Y-like, with cos(/?) =0.1. 


B. Benchmark scenarios 

Taking in account the scan of the parameters space, we chose the benchmark scenarios, where their main features 
can be recapped as follows: 

• Scenario lb: 2HDM-III as 2HDM-I, with the couplings <pff given by P^hdm-iii = #2HDM-i + A g and cos(/3 — 

«) = 0.5, x u kk = 1.5 (k= 2,3), X %2 = 1-8, x£j = l- 2 , X 23 = 0-2, V 22 = 0.5, xls = l- 2 , x | 3 = 0.1, m A = 100 GeV 
and m H ± =110 GeV, taking X and Y located in the blue region of the left panel from Figure [U 

• Scenario Ila: 2HDM-III as 2HDM-II, namely, the couplings <pff given by ^hdm-iii = ^ 2 HDM-ii + Ag and 
cos(/3 - a) = 0.1, X22 = 0.5, X33 = 1-4, X ^ 2 = 2 , X d 33 = 1-3, X ^ 3 = “0.53, x ^ 3 = 0.2, x | 2 = 0.4, x £ 33 = l- 2 , 
X23 = 0.1, rrtA = 100 GeV and ra#± = 110 GeV, taking X and Y allowed in the right panel of the Figure [lj 


6 Current low energy constraints on the Higgs boson masses have been studied very recently | 57 l, |58] . 
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FIG. 1. The allowed region in the plane I vs 7, using the constraint Eq. m, which is obtained from the radiative inclusive 
decay B —>> X s j. We obtain the Scenario lb, which is shown in the left panel, with 0.1 < cos(/3 — at) < 0.5, \kk = 1-5 (k= 2,3), 
X22 — 1.8, X33 — 1-2, X 23 * ~ 0.2, X22 — 0-5, X33 — 1-2, X23 =0.1. For Scenario Ha and Y, the allowed region is given in the 
right panel with cos(/3 — a) = 0.1, X22 — 0-5, X33 — 1-4, X22 — 2, X33 — 1-3, X23 — —0.53, X23 — 0-2, X22 — 0-4, X33 — 1-2, 
X23 = 0.1. For both cases m h = 125 GeV, 130 GeV < m H < 300 GeV, 100 GeV < m A < 250 GeV, 110 GeV < m H ± < 200 
GeV. 


• Scenario Y: 2HDM-III as 2HDM-Y, namely, the couplings <pff given by ^hdm-iii = #2HDM-y + ^9 an d 


cos(/3 - a) = 0.1, X 22 = 0.5, xh = 1.4, X d 22 = 2, X 33 = 1-3, xh = “0.53, xtz = 0.2, X £ 22 = 0.4, X | 3 = 1-1, 
X 23 = 0.1, rriA = 100 GeV and rn H ± = 110 GeV, taking the same X and Y for the Scenario Ha. 


Hereinafter, we only simulated benchmarks where cr.BR^cf) —» bs) (cross section of the charged current production 
v e (j)q multiplied by Branching Ratio of the channel decay 0 —» bs+ c.c., with (0 = ft, if)), are more than 0.15 fb 
so that, for an integrated luminosity of 100 fb -1 , we can start with at least 15 events. Finally, when producing 
differential spectra of physical observables, we will concentrate on three 2HDM scenarios where the number of Higgs 
signal events in the bs + c.c. mode are large enough in order to be able to appreciate the underlying dynamics. 


IV. NUMERICAL ANALYSIS 

In this section, we describe first the production of Higgs signal. We then discuss the most important SM backgrounds 
and the different kinematics selections on the simulated events. 


A. Higgs bosons signals 

We consider the leading production processed of Higgs boson: z/ e 0g, where (j) = ft and H while q is a light-flavor 
quark (i.e., u,d,s,c). We assume that 0 is dominantly decaying into bs (plus charge conjugation). So both of our 
signals, the lighter Higgs as well as the heavier Higgs one, contain three jet (one is forward and two are central), 
missing (transverse) energy and no-lepton. Out of the two central jets, one is 6-tagged and the other is a light-flavor 
jet. We estimated the parton level signal cross sections with flavor-violation within the 2HDM-III by using CalcHEP 
|59j . This implementation also calculates the BRs of the Higgs boson 0 into bs. For estimating the cross sections at 
the LHeC |44 l60l - [64{ . we consider an electron beam, of energy E e ~— 60 GeV and a proton beam of energy E p = 7000 
GeV, corresponding to a center-of-mass energy of approximately y/s = 1.296 TeV. The integrated luminosity is 100 
fb -1 . To estimate the event rates at parton level we applied the following basic pre-selections: 

p q T > 15 GeV, A R(q, q) > 0.4 (14) 


7 The charged-current production is approximately 5 times larger than the neutral current production. Moreover, the neutral current 
production contains an electron and, since we are vetoing leptons in this particular analysis, we consider only the charged processes. 
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with A R = A 77 2 + A</> 2 , where 77 and p are the pseudo-rapidity and azimuthal angle, respectively. We take m t = 173.3 
GeV as the top-quark pole mass. We set the renormalization and factorization scale at the Z-boson mass (which is 
approximately the momentum transfer scale for the signal) and adopt CTEQ6L [65j as Parton Distribution Functions 
(PDFs), with a s (the strong coupling constant) evaluated consistently at all stages (PDFs, hard scattering and decays). 

Considering the latter, we calculate in the allowed regions given above in the Figure [lj the event rates ( a.BR.L ) 
at parton level for the neutral Higgs bosons h and H in the Scenarios lb, Ha and Y, respectively, considering both 
luminosities of 100 fb -1 (left panel) and 1000 fb -1 (right panel), which are shown in the Figures [2H71 One can see 
that the blue region contains the best benchmark points for all scenarios. We show that the most optimistic is in 
fact Scenario lb for both Higgs bosons h and 77, which reach events rates of order 500-1300 (5000-13000) with an 
integrated luminosity of 100 fb -1 (1000 fb -1 ), although Scenarios Ha and Y also have some interesting benchmark 
points where one can obtain 17 events rates at the same luminosity. The Tab. un shows the benchmark points that 
we select as interesting for studies at the LHeC. There are twenty-seven in total, obtained by taking the same three 
different values of the 77 mass (uih = 130, 150, 170 GeV) in correspondence to nine different configurations of the 
other parameters.The product of cross sections times the relevant BRs (a.bs) are shown in Tab. HT1 




a.BR.L 

5000 
4000 
3000 
2000 
1000 


X 


X 


FIG. 2. Event rates ( cr.BR.L ) at parton level for the neutral Higgs boson h , where L is the integrated luminosity. We show 
Scenario lb for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider rrih — 125 GeV. 



a.BR.L 

400 

350 

300 

1 250 
200 
150 
100 
50 


FIG. 3. Event rates ( a.BR.L ) at parton level for the neutral Higgs boson 17, where L is the integrated luminosity. We show 
Scenario lb for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider uih = 130. 


Restricting ourselves to the points for which the inclusive event rates are most optimistic, all estimated by taking 
all the light-flavor quarks, the 5-quark and the gluon as fluxes inside the proton and upon considering appropriate 
flavor-mixing where appropriate, we have then proceeded as follows. The top-quark and W-boson were allowed to 
decay freely as implemented in PYTHIA [66j. Following this, it was recognized that the signal processes have unique 
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FIG. 4. Event rates ( a.BR.L ) at parton level for the neutral Higgs boson h , where L is the integrated luminosity. We show 
Scenario Ila for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider mu = 125 GeV. 
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FIG. 5. Event rates ( a.BR.L ) at parton level for the neutral Higgs boson H : where L is the integrated luminosity. We show 
Scenario Ila for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider ttih = 130. 



FIG. 6. Event rates ( a.BR.L ) at parton level for the neutral Higgs boson h : where L is the integrated luminosity. We show 
Scenario Y for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider mu — 125 GeV. 


kinematic profiles. In particular, the final state quark transverse momentum is less than the mass of the vector bosons, 
its energy is very high with a small angle with respect to the beam directions (i.e., high forward rapidity). This will 
serve as guidance in our event selection. However, before proceeding further in this direction, we have to acknowledge 





















9 



FIG. 7. Event rates ( a.BR.L ) at parton level for the neutral Higgs boson if, where L is the integrated luminosity. We show 
Scenario Y for 100 fb -1 (left panel) and 1000 fb -1 (right panel). We consider tyih — 130. 
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X 

Y 
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bs 

125 GeV 
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bs 

= 130 GeV 

a.bs 

m H 

bs 

= 150 GeV 

a.bs 

m H 

bs 

= 170 GeV 

a.bs 

Ib35 

28 

10 

28 

15.66 

6.392 

51.8 

1.209 

51.6 

0.30 

1.58 

0.117 

Ib47 

30 

5 

30 

16.14 

3.086 

48.2 

10.983 

48.0 

0.127 

1.80 

0.839 

Ib57 

44 

5 

44 

17.58 

11.861 

38.6 

5.14 

38.4 

2.303 

3.68 

0.137 

Hall 

20 

2 

20 

1.42 

1.055 

25.2 

0.097 

25.0 

0.091 

24.8 

0.085 

IIal4 

26 

2 

26 

1.44 

1.651 

26.0 

0.059 

25.8 

0.054 

25.6 

0.049 

IIa26 

36 

1 

36 

1.46 

1.621 

26.4 

0.045 

26.2 

0.042 

26.0 

0.038 

Yall 

20 

2 

-2 

1.42 

1.084 

25.2 

0.062 

25.0 

0.059 

24.8 

0.054 

Yal2 

22 

2 

-2 

1.44 

1.078 

25.6 

0.057 

25.4 

0.053 

25.2 

0.048 

Yal4 

26 

2 

-2 

1.46 

1.441 

26.0 

0.057 

25.8 

0.053 

25.6 

0.049 


TABLE II. Parameters for few optimistic benchmark points in the 2HDM-III as a 2HDM-I, -II and -Y configuration. Here bs 
stands for BR(^ —»> bs + bs), in units of 10 — 2 , where <f> — f, if, while a.bs stands for the cross section multiplied by the above 
BR as obtained at the LHeC in units of fb. We have analyzed only the benchmarks where the a.bs is greater than 0.15 fb, so 
that at least 15 events are produced for 100 fb -1 . 


at this point that these processes and their kinematic features to discover generic Higgs bosons have been studied 
since long IsizLIsi. Further parton level studies have been performed within the SM recently in [69j. In the context 
of BSM physics, cross section estimates while taking into account Next-to-Leading Order (NLO) correction factors 
have been performed in [70] but no signal and background studies have been reported therein. In fact, dedicated 
simulations at the event generator level have not been done extensively and we focus here on these aspects, most 
importantly, with the intent of detecting two Higgs bosons simultaneously in novel flavor-violating modes. 


B. Backgrounds 

There are mainly two groups of SM noise to our Higgs signals. The charged-current backgrounds, utb, vbbj , vb2j, 
vZj, and the photo-production ones, e~bbj , e~tt. For estimating the cross sections of these SM backgrounds, we have 
used the same pre-selections like for the signal, Eq. m, and identical conventions and parameter sets. The expected 
number of events for 100 fb -1 of integrated luminosity are given in the third column of Tab. IIII1 


C. Signal-to-background analysis 

We passed the CalcHEP v3.4.7 [59j generated parton level event on to PYTHIA v.6.408 (66[, which handles the parton 
shower (both initial and final), hadronization, heavy hadron decays etc. The final state radiation smears the four- 
momentum of the jets, thus the invariant mass of the Higgs boson signal is less than the actual values considered in 
the event. We also took the experimental resolutions of the jet angles and energy using the toy calorimeter PYCELL, 
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Proc 

SimEvt 

RawEvt 

a 

b 

c 

d 

e 

f 

g 

h 

i 

5 

Ib35 

100K 

639.2 

447.6 

177.3 

117.1 

97.4 

93.8 

37.8 

31.7 

25.4 

15.8 

1.2(3.8) 

Ib47 

100K 

308.6 

216.8 

85.1 

56.2 

47.1 

45.5 

18.4 

15.6 

13.0 

8.1 

0.62(2.0) 

Ib57 

100K 

1186.1 

833.7 

325.7 

215.5 

180.6 

173.9 

70.3 

59.1 

49.3 

31.1 

2.4(7.5) 

Hall 

100K 

105.5 

74.3 

29.1 

19.2 

16.0 

15.4 

6.3 

5.3 

4.4 

2.8 

0.21(0.70) 

IIal4 

100K 

165.1 

116.1 

45.2 

30.0 

25.4 

24.4 

9.7 

8.3 

6.9 

4.4 

0.33(1.05) 

IIa26 

100K 

162.1 

114.4 

44.7 

29.5 

24.5 

23.6 

9.5 

8.1 

6.8 

4.3 

0.33(1.03) 

Yall 

100K 

108.4 

76.3 

29.8 

19.6 

16.4 

15.8 

6.4 

5.4 

4.6 

2.9 

0.22(0.70) 

Yal2 

100K 

107.8 

76.2 

29.6 

19.5 

16.3 

15.7 

6.3 

5.4 

4.5 

2.8 

0.21(0.67) 

Yal4 

100K 

144.1 

101.7 

39.8 

26.0 

21.7 

20.8 

8.2 

7.0 

5.9 

3.8 

0.29(0.92) 

utb 

100K 

50712.1 

28338.4 

15293.7 

9845.0 

8144.2 

7532.7 

2982.1 

2058.0 

652.2 

139.6 


ubbj 

560K 

14104.6 

6122.8 

3656.7 

1858.5 

1787.1 

1650.1 

257.5 

152.5 

85.2 

15.1 


ub2j 

90K 

18043.1 

8389.2 

3013.0 

1691.5 

1445.5 

1373.7 

389.5 

206.1 

77.2 

11.3 

.8=170.8 

u3j 

300K 

948064.2 

410393.4 

15560.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

>8=13.1 

ebbj 

115K 

256730.1 

55099.8 

36353.6 

12659.8 

1432.0 

200.7 

54.1 

24.8 

18.0 

4.5 


ett 

130K 

783.3 

685.0 

384.5 

265.9 

179.3 

26.2 

11.6 

10.5 

3.9 

0.3 



TABLE III. Expected number of events after different combinations of cuts for signal and backgrounds at the LHeC with 100 
fb -1 integrated luminosity for m^=125 GeV. SimEvt stands for the actual number of events analyzed in the Monte Carlo 
simulations. RawEvt stands for the number of events with only the generator-level cuts imposed; for the signal as well 
as for background, these are calculated from the total cross section times branching ratio. In the final column we mention the 
significances(<S) defined as S = S/y/B, where signal events S, background events B for 100 fb -1 of data after all cuts mentioned 
in the “i” column. The number in the parenthesis in the final column represent the significances for 1000 fb -1 . 


in accordance with the LHeC detector parameters, given in PYTHIA. This has some non-trivial effect since we used 
the invariant mass to isolate the Higgs signal. In our study we considered the LHC type calorimeter for the LHeC. 
Although in reality this is not the case, for example, unlike ATLAS and CMS, the electro-magnetic and the hadronic 
calorimeter at the LHeC is not symmetric. However, since we are not doing detector simulation and also we are not 
considering cracks in the detectors, we applied symmetric large rapidity coverage for jets and leptons in our analysis. 
We expect that these assumptions hardly alter our findings. The detector parameters in the toy calorimeter module 
PYCELL are set according to the LHeC detector [62]. Specifically, we assume large calorimeter coverage \rj\ < 5.5, with 
segmentation (the number of division in 77 and </> are 320 and 200 respectively) A 77 x Ac/) = 0.0359 x 0.0314. Further, 
we have used Gaussian energy resolution [60] for both leptons (£ = e,/i) and jets (labelled as j), with 


A E 
~E~ 



(15) 


where a = 0.32, b = 0.086 for jets and a = 0.085, b = 0.003 for leptons and ® means addition in quadrature. We 
have used a cone algorithm for the jet-finding algorithm, with jet radius A R(j) = yj A 77 2 + A (j) 2 = 0.5. Calorimeter 
cells with >5.0 GeV are considered to be potential candidates for jet initiators. All cells with >1.0 

GeV were treated as part of the would-be jet. A jet is required to have minimum summed Erp min > 15 GeV and 
the jets are ordered in Et . Leptons {£ = e, fj) are selected if they satisfy the requirements: > 15 GeV and 

I?/1 < 3.0. In our jet finding algorithm we include leptons as parts of jets. Finally we separate them, putting some 
isolation criterion as follows: if we find a jet near a lepton, with AR(j — £) < 0.5 and 0.8 < E° t /Ej, < 1.2, i.e. if 
the jet Et is nearly identical to that of this lepton, the jet is removed from the list of jets and treated as a lepton. 
However, if we find a jet within A R(j — £) < 0.5 of a lepton, whose Et differs significantly from that of the lepton, 
the lepton is removed from the list of leptons. This isolation criterion mostly remove leptons from b or c decays. We 
reconstructed the missing (transverse) energy (pr) from all observed particles and shown in left panel of Fig. [9] We 
have also calculated the same from the energy deposition in the calorimeter cells and found consistency between these 
two methods. Only jets with \rf\ < 2.5 and E° T > 15 GeV “matched” with a b —flavored hadron (B— hadron), i.e. 
with A R(j,B — hadron) < 0.2 is considered to be “taggable”. We assume that these jets are actually tagged with 
probability = 0.50. We also adopted mis-tagging of non —b jets as b —jets and treated c— jets differently from the 
gluon and light-flavor jets. A jet with \ < 2.5 and E° T > 15 GeV matched with a c— flavored hadron (C— hadron, 
e.g., a D —meson or A c —baryon), i.e., with A R(j,C — hadron) < 0.2, is again considered to be taggable, with (mis- 
)tagging probability e c = 0.10. Jets that are associated with a r—lepton, with A R(j,r) < 0.2, and all jets with 
|t 7 j | > 2.5, are taken to have vanishing tagging probability. All other jets with E° T > 15 GeV and |? 7 J | < 2.5 are 
assumed to be (mis-)tagged with probability e Uj d,s,g = 0.01. These efficiencies follow recent LHeC analyses f69|. 

The analysis strategy has been adopted from earlier work of some of us [71 ]. In particular, we have exploited a 
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SimEvt 

RawEvt 

A 

B 

C 

D 

E 

F 

G 

H 

I 

5 

Ib35 

100K 

120.9 

87.1 

34.1 

26.9 

22.5 

21.6 

7.5 

6.1 

5.3 

3.4 

0.28(0.88) 

Ib47 

100K 

1098.3 

790.3 

307.1 

243.9 

204.6 

195.7 

68.5 

56.1 

48.6 

31.3 

2.6(8.1) 

Ib57 

100K 

514.0 

371.2 

144.8 

115.0 

96.0 

92.0 

31.7 

25.8 

22.7 

14.3 

1.2(3.7) 

Hall 

100K 

9.7 

6.8 

2.7 

2.1 

1.8 

1.7 

0.6 

0.4 

0.3 

0.2 

0.02(0.05) 

IIal4 

100K 

5.9 

4.2 

1.7 

1.3 

1.1 

1.0 

0.4 

0.3 

0.2 

0.1 

0.01(0.02) 

IIa26 

100K 

4.5 

3.1 

1.3 

1.0 

0.8 

0.8 

0.3 

0.2 

0.1 

0.1 

0.01(0.02) 

Yall 

100K 

6.2 

4.4 

1.8 

1.4 

1.1 

1.1 

0.4 

0.3 

0.2 

0.1 

0.01(0.02) 

Yal2 

100K 

5.7 

4.0 

1.6 

1.3 

1.0 

1.0 

0.3 

0.2 

0.2 

0.1 

0.01(0.02) 

Yal4 

100K 

5.7 

4.0 

1.6 

1.3 

1.0 

1.0 

0.3 

0.2 

0.2 

0.1 

0.01(0.02) 

utb 

100K 

50712.1 

28338.4 

15293.7 

10976.4 

9092.4 

8393.6 

2550.9 

1565.5 

617.9 

113.7 


vbbj 

560K 

14104.6 

6122.8 

3656.7 

2145.5 

2062.1 

1902.9 

266.6 

141.0 

87.5 

14.4 


ub2j 

90K 

18043.1 

8389.2 

3013.0 

2053.6 

1734.0 

1650.1 

402.8 

143.7 

64.5 

8.1 

.8=147.8 

u3j 

300K 

948064.2 

410393.4 

15560.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

v / 8=12.2 

ebbj 

115K 

256730.1 

55099.8 

36353.6 

16838.4 

1826.6 

284.1 

56.4 

31.6 

22.6 

11.3 


ett 

130K 

783.3 

685.0 

384.5 

280.8 

190.8 

27.8 

10.9 

9.3 

3.9 

0.3 



TABLE IV. Same as Tab. |m] but for rriH= 130 GeV. The criterion for jets and 5-tagging are the same, so that the number of 
events in column A and B are the same for all SM backgrounds. 


Proc 

SimEvt 

RawEvt 

A 

B 

C 

D 

E 

F 

G 

H 

I 

<S 

Ib35 

100K 

30.0 

23.3 

9.1 

8.2 

6.9 

6.5 

1.5 

1.3 

1.2 

0.8 

0.10(0.33) 

Ib47 

100K 

12.7 

9.9 

3.8 

3.4 

2.9 

2.7 

0.6 

0.5 

0.5 

0.3 

0.04(0.12) 

Ib57 

100K 

230.3 

179.6 

69.3 

62.6 

52.6 

49.9 

11.7 

10.1 

9.1 

6.4 

0.83(2.62) 

Hall 

100K 

9.1 

6.9 

2.7 

2.4 

2.0 

1.9 

0.4 

0.4 

0.3 

0.2 

0.026(0.08) 

IIal4 

100K 

5.4 

4.1 

1.6 

1.4 

1.2 

1.1 

0.3 

0.2 

0.2 

0.1 

0.013(0.04) 

IIa26 

100K 

4.2 

3.2 

1.3 

1.1 

0.9 

0.9 

0.2 

0.1 

0.1 

0.1 

0.013(0.04) 

Yall 

100K 

5.9 

4.5 

1.8 

1.6 

1.3 

1.2 

0.3 

0.2 

0.2 

0.1 

0.013(0.04) 

Yal2 

100K 

5.3 

4.0 

1.6 

1.4 

1.2 

1.1 

0.3 

0.2 

0.2 

0.1 

0.013(0.04) 

Yal4 

100K 

5.3 

4.0 

1.6 

1.4 

1.2 

1.1 

0.3 

0.2 

0.2 

0.1 

0.013(0.04) 

istb 

100K 

50712.1 

28338.4 

15293.7 

11810.9 

9808.7 

9039.0 

751.7 

476.8 

194.5 

32.3 


vbbj 

560K 

14104.6 

6122.8 

3656.7 

2395.6 

2300.1 

2120.8 

199.3 

112.4 

70.8 

12.4 


vb2j 

90K 

18043.1 

8389.2 

3013.0 

2427.2 

2030.3 

1933.1 

234.2 

83.7 

41.0 

6.3 

B=60.1 

v3j 

300K 

948064.2 

410393.4 

15560.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

VB= 7.7 

ebbj 

115K 

256730.1 

55099.8 

36353.6 

21280.9 

2270.8 

385.6 

36.1 

24.8 

20.3 

9.0 


ett 

130K 

783.3 

685.0 

384.5 

291.5 

199.0 

29.1 

3.5 

3.0 

1.2 

0.1 



TABLE V. Same as Tab. [TV|but for m H = 150 GeV. 


simple cut-based method for signal enhancement and background rejection. We have chosen the following selections 
and applied them cumulatively for the signal from h (H). 

• a(A): We first selected events containing at least three jets (same). The distribution of the number of jet ( Nj et ) 
is shown in the left panel of Fig. [51 For the lighter Higgs, ft, in all the signal benchmarks the efficiencies are 
approximately 70%(For heavier Higgs boson iL, with mass of 130, 150 and 170 GeV the efficiencies are 71%, 
76% 81% respectively, as two of the jets are directly coming from the corresponding heavier Higgs boson and 
hence the efficiencies are larger.). Further, tb has efficiency of approximately 56%, 2bj and 3 j events display 
approximately 45% whereas the 52 j one has approximately 47%. Efficiencies are generally lower for the photo¬ 
production channels: e2bj is approximately 20% (the sharp fall is due to the isolation criterion) whereas for e2t 
the jet efficiencies are higher and due to the presence of two 5-quarks leads to two 5-quarks and the probability 
of having two jets from W-boson decay itself is approximately 91%, thus, out of 4-jets in 91% events, the 
probability of having at least three energetic jets is reduced by 4%, which leads to approximately 87% efficiency. 
However, as we will see, the presence of the electron in the photo-production modes leads to the corresponding 
backgrounds to be very low. 


8 Unless mentioned otherwise, all the efficiencies quoted hereafter are given with respect to the previous selection. 
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RawEvt 
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I 

5 

Ib35 

100K 

11.7 

9.6 

3.7 

3.5 

3.0 

2.8 

0.5 

0.4 

0.4 

0.3 

0.053(0.17) 

Ib47 

100K 

83.9 

69.2 

26.7 

25.5 

21.5 

20.2 

3.6 

3.1 

3.0 

2.2 

0.39(1.23) 

Ib57 

100K 

13.7 

11.2 

4.3 

4.1 

3.4 

3.2 

0.6 

0.5 

0.5 

0.4 

0.07(0.22) 

Hall 

100K 

8.5 

7.0 

2.7 

2.5 

2.1 

2.0 

0.3 

0.3 

0.3 

0.2 

0.035(0.11) 

IIal4 

100K 

4.9 

4.1 

1.6 

1.5 

1.3 

1.2 

0.3 

0.17 

0.16 

0.12 

0.021(0.07) 

IIa26 

100K 

3.8 

3.1 

1.2 

1.1 

0.9 

0.9 

0.1 

0.1 

0.1 

0.1 

0.02(0.06) 

Yall 

100K 

5.4 

4.4 

1.7 

1.6 

1.4 

1.3 

0.2 

0.2 

0.2 

0.1 

0.02(0.06) 

Yal2 

100K 

4.8 

4.0 

1.5 

1.4 

1.2 

1.1 

0.2 

0.2 

0.2 

0.1 

0.02(0.06) 

Yal4 

100K 

4.9 

4.0 

1.6 

1.5 

1.2 

1.1 

0.2 

0.2 

0.1 

0.1 

0.02(0.06) 

utb 

100K 

50712.1 

28338.4 

15293.7 

12381.7 

10299.7 

9465.2 

209.7 

144.5 

75.9 

13.2 


ubbj 

560K 

14104.6 

6122.8 

3656.7 

2568.2 

2465.8 

2272.4 

103.7 

60.8 

37.4 

8.7 


vb2j 

90K 

18043.1 

8389.2 

3013.0 

2744.8 

2278.1 

2171.4 

99.5 

40.0 

25.2 

5.3 

B= 31.7 

u3j 

300K 

948064.2 

410393.4 

15560.9 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

YB= 5.6 

ebbj 

115K 

256730.1 

55099.8 

36353.6 

25010.7 

2638.4 

453.3 

29.3 

18.0 

11.3 

4.5 


ett 

130K 

783.3 

685.0 

384.5 

298.8 

204.5 

29.9 

1.0 

0.8 

0.5 

0.0 



TABLE VI. Same as Tab. [IV|but for rriH= 170 GeV. 




FIG. 8. Scenario lb with the parameter X = 28 = Z and Y = 10. The number of jets ( Nj e t ) in the left panel and the 
number of fr-tagged jet with the inclusion of mis-tagging (Nb-tag) in the right panel for signal (rriH= 150 GeV) and all the SM 
backgrounds. For other signal events, the distribution profiles are very similar, except the fact that the number of jets as well 
as 6-tagged jets is slightly larger for heavy Higgs bosons. See the fourth(hfth) column for their actual efficiencies with selections 
applied in a(b) and A(B) for lighter (heavier) Higgs bosons, respectively. 


. b(B): We demanded at least one 6-tagged jet with the inclusion of proper mis-tagging (same). The distributions 
of the number of 6-tagged jets (Nb-tag) are shown in the right panel of Fig. [HI For the lighter Higgs, ft, in all the 
signal benchmarks the 6-tagging efficiencies are approximately 40%. In fact, all our signal benchmarks contain 
at least one 6-quark and, since we adopted 65=0.50, the 10% lowering is quite realistic and due to the fact 
that not all 6-quarks in the signal are eligible for the 6-taggable criterion adopted in our analysis. For heavier 
Higgs signals, benchmarks show similar efficiencies and changes are less than 1% for three masses. In case of 
66, the events containing at least one 6-tagged jet are approximately 54%, a rate greater than the signal ones, 
as these background topologies contain at least two 6-quarks at the parton level so that, due to combinatorics 
(other than mis-tagging a light-flavor quark-jet from LF-boson decays), the probability of one 6-tag is more. The 
probability of 2 bj is approximately 60%, approximately 6% larger than 66, due to, unlike 66, the presence of one 
6-quark and one light-flavor jet in the hard processes. Further, b2j efficiencies are similar to or little less than 
those of the signals due to the fact that the taggable rapidity is more central, where the jets are more likely to 
be forward in the basic hard processes. The efficiency of e26 is approximately 2% larger than for 66 mainly due 
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to the contributions from mis-tagging one extra hadronic W, i.e, two extra jets. Finally, the efficiency of 3 j is 
approximately 4%, which seems consistent with the expected mis-tagging rates (1.0% for light-flavor jets and 
10% for c-jets) in presence of combinatorics. 

• c(C): We demanded at least two central jets, with pp > 30 GeV (pr > 25, 20 and 15 GeV for mn— 130, 150 
and 170 GeV, respectivel}@) in the pseudo-rapidity range |^| < 2.5. One of the central jet must be a 6-tagged 
jet and we demanded only one 6-tagged jet (same). For the lighter Higgs, ft, all the signal benchmarks as well 
as the tb noise survive at the rate of approximately 65% since all these processes naturally have three jet in 
their events. Further, e2bj and 2 bj are reduced by approximately 35% and 50%, respectively, mainly due to 
demanding, with respect to the b(B) case above, of one 6-tagged jet only. The diagrams of tb and e2 1 reveal 
that these two backgrounds are more central, because of the presence of one t-quark in the central region. In 
fact, the efficiency is larger in e2t and is mainly due to the contributions from the additional t-quark. Although 
in 2 bj the probability in presence of one 6-tagged in the central region is large the overall efficiency is reduced 
to 12% due to (partly) the possibility of more than one 6-tagged jets whereas for Wb the value is mainly due 
to the hadronic branching fraction and also that the 5-tagged jet is not necessary central. In case of e2bj the 
efficiency is 35%: this noise suffers mainly due to the centrality criterion. For 3j, none of the events survived 
this selection criterion. The efficiencies pattern discussed above are similar for the heavier Higgs boson, ift, yet 
recall that here we have used slightly softer selections on the transverse momentum. Thus, the efficiencies are 
increasing with a softer pr selection for both signals and backgrounds. 

. d(D): The missing transverse energy cut pr > 20 GeV is first applied (same). The relevant distribution is 
shown in the left panel of Fig. [9j For all the signal benchmarks (lighter as well heavier Higgs bosons), tb and 52 j 




FIG. 9. The missing energy (pr) (left panel) and rapidity (rjj f ) (right panel) profile of the forward jet for signals and SM 
backgrounds. The pr distributions for all other signal benchmarks as well as the tb noise are not shown as they are very 
similar to the signal distributions of ttih— 150 GeV for Scenario lb with X — Z — 28 and Y — 10 (shown in thick solid), 
whereas the thin solid is for rrih= 125 GeV for Scenario la with X — Z — 28 with Y — 10. The rapidity distributions profile 
for 771^ = 130(170) GeV is very close to the 777^=125 GeV(777i7 = 130 GeV) case shown in thin solid, except that for massive 
Higgs the peaks shift towards the left. Also the corresponding rapidity distribution profile for e2bj is somewhat similar to the 
777 /i= 125 GeV signal case. 

the efficiencies are approximately 83%. The sources of neutrinos and the event structures of these two processes 
are very similar, except for the fact that the top-quark decays produce either neutrinos (which are then pure 
sources of missing energy but in such a case, owing to the selection a(A), they are largely removed) or quarks 
(where the smearing of jets and track mis-measurements are the main sources of missing energy). For 2 bj the 
efficiencies are approximately 96%. This selection is crucial to suppress the photo-production processes: e2bj 
and e2t. In case of e2bj only 12% of the events survive in fact. For e2t the presence of two VF-bosons and their 


Since the cross sections become smaller with increasing m#, we lowered the central jet pr cuts. 
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decays into leptonic modes would be the sources of missing energy, so that this noise is not affected very much 
by our constraints. 

. e(E): A lepton (e or (i) veto for pr > 20 GeV and 77 < 3.0 is applied here (same). For the lighter Higgs signal 
benchmarks, the efficiency for this selection is approximately 96%, as only 4% of the events contain at least one 
lepton which is coming from the semileptonic decays of H-hadrons and D-mesons and which passes the isolation 
criterion above. The efficiencies for 2bj and b2j are rather close to the signal benchmarks, due to these processes 
also not having prompt leptons in their events. The tb channel has an efficiency of 93%, hence approximately 
3% less than the signal, as here the top-quark decays can lead to one bottom-quark and, if the hard-processes 
bottom quarks are more central, the requirement in c(C) is satisfied and there is no problem in having a central 
and high-pT lepton from a VF-boson decay. In the photo-production processes, e2 bj and e2t both contain hard 
leptons, so only 15% of the events survive our lepton veto. For heavier Higgs masses, the efficiencies are quite 
similar to the lighter Higgs boson mass case. The only difference is that, since the applied pr threshold is lower 
for heavier Higgs masses, the probability of having a lepton in the event is higher, thus the veto efficiencies are 
somewhat smaller, except for the two photo-production processes. 

• f(F): In the central region, defined above via c(C), we reconstruct the invariant mass of one 6 -tagged jet with any 
of the other jets, M^. Amongst these, we have chosen the best combination, i.e., where, the absolute difference 
|Mbj — M h ( H ^ \ is minimized. We call this di-jet combination the candidate light (heavy) Higgs boson signal. In 
order to select the latter, we have kept events within a 15 GeV mass window centered around the corresponding 
Higgs boson masses. The distributions of M^j are shown in the left panel of Fig. [lOl It seems that the di-jet 
invariant masses of the BP-Ia30 signal benchmark (with 777^=125 GeV) has its peak around 115 GeVF^l. The 
distribution for tb also has a peak around that of the 777^=125 GeV signal. However, the combinatorics is 
significant and this shows in their efficiencies, which are approximately 40% for both. The distribution of b2j 
is flat as there is no correlations for the correct di-jet candidates. Also note that b2j has TF-boson exchange 
resonant diagrams, so a probability in principle exists for a di-jet invariant mass peak at Mw-, however, this 
is very small, mainly due to low mis-tagging efficiencies and the centrality criterion. Further, also in case of 
2bj, where the Z-boson is present resonantly in the diagrams, the Z bb decay combined with high tagging 
efficiencies allows for the appearance of a di-jet peak at 80 GeV (approximately 10 GeV less than Mz due to jet 
energy smearing): see the left-panel of Fig. [TUI In case of e 2 £, like tb, one has also correlated di-jet candidates, 
but the energy scale is higher, so the peak is shifted to higher masses. The efficiency is approximately 45%, a 
little larger than the signal and tb ones. The distributions of b2j and e2bj are flat and the efficiencies are the 
same, approximately 28%. This particular selection suppresses 2 bj events more severely though, at a rate which 
is approximately 15%. For heavier Higgs bosons the distributions show rapid falls and so, by applying the mass 
window cuts, only the left part of the distributions contributes. This shows in their signal efficiencies, which are 
approximately 34%, 23% and 18% for ttt,#= 130, 150 and 170 GeV, respectively. The SM backgrounds do not 
show up in distributions at large invariant masses, thus for heavier Higgs mass combined with the same mass 
window selection suppresses more the backgrounds. As an example, in case of e 2 £, which produces somewhat 
higher invariant masses than all other SM backgrounds, the efficiencies drops to 40%, 12 % and 4% for ttt,#= 130, 
150 and 170 GeV, respectively. In case of tb the efficiencies (see Tabs. ITTT1 HVl lyland lVTj) drop to from 30%, 
8 % and 2%, respectively. In case of 2 bj the values are 14%, 10% and 5%. For e2bj, one has 20%, 10% and 7%, 
respectively. Finally, for b2j, these are 25%, 12 % and 5%, respectively. 

• g(G): We demanded the remaining leading jet in the event to have pr > 25 GeV, with —5.5 < 77 < —0.5 (—1.0) 
(these values are chosen by seeing the distribution, see the right panel of Fig. [9]) and termed it as the forward 
jet (jf). This forward jet lies very close to the direction of the incoming proton, like in Vector Boson Fusion 
(VBF) processes for Higgs production. In contrast to VBF though, instead of a jet with large rapidity gap with 
respect the the forward jet (a backward jet then), in our signal we have a neutrino. The more massive the Higgs 
is, the less energy remains for the forward jet so as to lay close to the proton direction, i.e., at larger rapidity. 
This reflects in the right panel of Fig. [9] The thick (thin) solid curve corresponds to 777 ,#=150 GeV (777^=125 
GeV). For a lighter Higgs boson, 777^=125 GeV, the efficiency is approximately 80%. For e2t the efficiency is 
almost 90%, twice that of e2bj, as there is more than twice a probability to have a forward-jet from top-quark 
decays. 

. h(H): The di-jet invariant mass of the Higgs boson candidates with the forward tagged jet, which is essentially 
the overall energy scale of the hard scattering, is asked to comply with the following requirements: rrihj f ( rriHj f ) 
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The peaks always show up to the left side of the actual masses due to jet energy smearing and the shift also depends on the jet-cone 
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FIG. 10. The di-jet invariant mass, made up by one 5-tagged and one light-flavor jet, producing Higgs candidates, = M^j 
(left panel) and the three-jet invariant mass, i.e., the previous two jets combined together with the forward jet, (right 

panel). The mass peaks of the Higgs signals (M#) correspond to rrih= 125 (thin black) for Scenario la, rriH= 150 (thick black) 
and 170 (thin black) for Scenario lb from left to right. All these are using the parameters X — Z — 28 and Y =10. The 
distribution for 77777=130 is not shown but it lies in between 777^=125 and 77777 = 150. Among all SM backgrounds, only 2bj 
shows a prominent peak from the Z-boson. Notice that M^j f represents the overall energy scale of the hard-scattering. 


> 190 GeV (190, 210 and 230 GeV for ra#=130, 150 and 170 GeV, respectively). The distributions are shown 
in the right panel of Fig. [lOl For a lighter Higgs boson, 777^=125 GeV, except a few case£3, the efficiency is 
approximately 82%. This forward jet should not be a 5-jet though. So, in 55 and e25, where a forward 5-tag jet 
is more probable, the efficiencies are lower, approximately 32% and 37%. It is clear from the right panel of Fig. 
EH that the three-jet invariant mass distributions of 52j peak around 140 GeV or so. The same for 55, e25, 25j, 
e2bj and the Higgs signal with 777^=125 GeV, which show somewhere around 180 GeV. So, for 777^=125 GeV, 
the efficiency is around 82%. For the heavier Higgs boson, with ran = 150 (170) GeV, the distributions are 
shown in the thick (thin) solid curve in the right panel of Fig. [TOland peak around 220 (260) GeV. The selection 
cuts for these two Higgs bosons are 210 and 230 GeV, respectively. With our selection, for these two heavy 
Higgs signals, one suppresses more SM backgrounds than in the case of the Higgs signal with 777,^=125 GeV and 
rriH = 130 GeV. For example, in the tb case, the most dominant background, for 777^=125 GeV, and tt 7,#=130, 
150 andl70 GeV, the events which survived are approximately 652, 618, 195 and 80, respectively. For other SM 
backgrounds a similar pattern follows, thus the overall SM backgrounds can be strongly reduced. However, this 
overall background rejection will not help alone to have larger significances as the signal rate itself is suppressed 
with heavier Higgs masses (at the production level). 

• i(I): Finally, we required only one light-flavor jet in the central region (same). This selection is called ‘central jet 
veto’ and has severe impact on all processes having more jets in the central rapidity region, other than the Higgs 
candidate jets. Recall that our Higgs signal candidate jets, selected in f(F) above, are central: this is true for not 
only the signal, but also the dominant SM background, tb. For a lighter SM Higgs, see Tab. Hill approximately 
35-40% of the events have a central jet other than Higgs candidate jets, thus only 60-65% of the events survive. 
For 55, vb2j and v2bj the efficiencies are 22%, 18% and 14%, respectively. Among all the SM backgrounds, e25 
has a larger number of jets (see the distributions in the left panel of Fig. thus the probability of having a 
central jet is more, so that this selection suppresses this background severely, approximately by 93% (for all the 
Higgs cases, see Tabs. Hill HVHVl and ED). 

After the cumulative selections from a-i, discussed above, we find that, for the SM Higgs boson with 777^=125, the 
final number of events is around 15-30 only for Scenario lb and for large values of the parameters X = Z = 28(44) 


ii 


Recall that the efficiencies are relative to the previous selection, one can estimate the individual efficiencies from the respective distri¬ 
butions. 
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and Y = 10(5) respectively. The total SM background rate is approximately 170. The charged-current backgrounds, 
vtb, vbbj and z/62j, are the dominant ones and only 3% of the total background comes from ett photo-production. 
These rates lead to a maximum significance of approximately 2.4 (7.5)a with 100 (1000)/6 _1 integrated luminosity 
for Scenario lb with X = Z = 44 and Y = 5. For Scenario lb with X = Z = 28 and Y = 10, the significance is 
approximately 1.2 (3.8)a for 100 (1000)/6 _1 integrated luminosity. The significances for Scenario Y and Scenario Ha 
are less than 1. Thus, one can expect that Scenario lb with large value of X = 44 = Z and Y = 5 may be accessible 
through the SM-like Higgs boson signal already detected at the LHC. 

We also searched for the second CP-even neutral Higgs boson of our 2HDM-III, with masses m#=130, 150 and 
170 GeV. After the cumulative selections from A-I, the maximum number of signal events for mu —130 GeV is 
approximately 15 (30) and only for Scenario lb with X = Z=44 (30) with Y = 5. The total SM background is 
approximately 150. So the maximum significance is approximately 2.6 (8.1)a for 100 (1000)/6 _1 integrated luminosity. 
For the case m#=150 GeV, the number of signal events is approximately 7 and the SM background reduces to 
approximately 60: this leads to a significance of approximately 1.0(2.6)a with 100 (1000)/6 _1 integrated luminosity. 
For rriH =170 GeV, the raw event rate is approximately 80 to start with and, at the end, the count is only £3. The 
total SM backgrounds is approximately 30, which leads to a significance approximately 0.4 (1.23)a. 

The LHeC will be operational for about ten years and expected to accumulate a total integrated luminosity of 1000 
fb -1 of data. So, at the end of the run, we expect the SM Higgs boson will have 7.5cr (3.8a) significance for Scenario lb 
with X = Z=44 (28) with Y = 5(10) . For Scenario Ha and Scenario Y, with A=26 and Y = 2, the final significances 
are approximately la. It seems that in all scenarios of the model, large X are favorable. For heavy Higgs masses with 
mn =130 GeV, for Scenario lb with X=44 and 30 with Y = 5, the maximal significances are approximately 8.1 and 
3.7a respectively. For ra#=150 GeV, Scenario lb with A=44 and Y = 5, the final significance is 2.62a. For mu —170 
GeV, in the Scenario lb with A=30 and Y = 5, the final number of signal events is approximately 2. The estimated 
significances is 1.23a. Thus, for high enough Higgs masses, one might invoke the aforementioned multivariate analyses 
to have larger significances. 


V. CONCLUSIONS 

After the discovery of a SM-like Higgs boson at LHC, one is well motivated to look for more such states, which 
necessarily appear in BSM scenarios. Among the experimental facilities where more Higgs bosons can be searched for, 
one should list an ep collider which may be built at CERN, known as the LHeC. In our analysis we have considered 
a 2HDM-III with a four-zero Yukawa texture in three configurations, wherein both the SM-like Higgs boson and the 
heavier version of it can be accessible at the foreseen LHeC energy. We assumed that both of these states are decaying 
via a flavor-violating mode (</> -A bs). After a parameter scan, we have selected a few model benchmarks where the 
products of cross sections and flavor-violating BRs are large enough to produce sufficient events in which to look for 
both signatures. We studied the three-jet and missing energy channel, 3 j + fir, from the charged-current production 
of v e (j)qf, where qf is a forward jet with large rapidity and the other two jets come from the flavor-violating decay 
(j) -A bs. We demanded one central jet to be 6-tagged with the inclusion of the proper mis-tagging from light-flavor 
and gluon jets. We considered the most dominant SM backgrounds: charged-currents, vtb , vbbj, vb2j and v3j, 
and photo-production, e~bbj and e~tt. We performed a full hadron-level Monte Carlo simulation using CalcHEP as 
matrix element calculator, PYTHIA as parton shower/hadronization event generator and its PYCELL toy calorimeter in 
accordance with the LHeC detector parameters. We carefully implemented 6-tagging, including mis-tagging of ojets 
or light-flavor or gluon jets. 

The signals under consideration do not have leptons, so we applied lepton vetos. The charged-current production 
has naturally missing energy due to the presence of neutrinos but no charged lepton. However, the photo-production 
processes have leptons in them but no direct missing energy (except the mis-measurements from jets and smearing), 
thus the missing energy selection together with the lepton veto suppressed the photo-production backgrounds to a 
very large extent. 

The kinematics of the particular signals considered is very interesting from the fact that the Higgs boson is pro¬ 
duced in the central rapidity region and its decay daughters, one 6-jet and one light-flavor jet, are also central. We 
reconstructed the invariant mass of this two jets and selected events only for masses within a 15 GeV window around 
the respective Higgs masses of the signal benchmarks. This selection reduces the SM backgrounds to a large extent 
and the invariant mass ensures the selection of flavor-violating decays. For massive Higgs bosons, although the signal 
events becomes low with the mass window selection, background suppression is more efficient. 
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Note that our selection cuts applied above are not optimized. An increase of the luminosity is an easy solution from a phenomenological 
perspective. However, adopting multivariate analysis techniques must also be a better discriminator of signal from backgrounds. 
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As a next step of our selection, we identified the most energetic light-flavor forward jet (by seeing the rapidity 
profiles) and calculated the invariant mass with that jet together with the flavor-violating Higgs candidates jets. 
These three-jet invariant masses essentially give the overall energy scale of the hard scattering. Again, the more 
massive Higgs boson helps to suppress more SM backgrounds, in particular utb and ett, but the signal becomes 
smaller too. 

At the end, the most important cut, we applied a central jet veto, i.e., we required one light-flavor central jet only. 
This suppresses SM backgrounds with large multiplicity, for example, vtb and ett. 

After all the selections, with 100 fb -1 of data, we found that the SM Higgs boson, ft, would be detectable within 
the 2HDM-III in the scenario called in this work lb withX = Z = 44 or 30 with Y = 5, with approximately l-2a. 
The heavier neutral Higgs boson, H, with masses 130 GeV (150 GeV), would have 2 (l)cr significances for large X 
and only for Scenario lb. 

The LHeC will be operational for around ten years and so it is expected to accumulate a total integrated luminosity 
of 1000 fb -1 of data. So, in all the cases mentioned above, the final significances will be enhanced. At the end of the 
run, the 2HDM-III Like-IIa (Like-Y) the SM Higgs will have la. For ra#=130 GeV, in the Scenario lb with X=44 
(30) and Y = 5, the maximal significances are approximately 3.7 (8.1)cr. The maximal significances for mu —150 
GeV is 2 .6a for lb with X=44 and Y = 5. For ra#=170 GeV the final number of signal events is approximately 2, 
probably too little to be detected. However, it should be noted that we have adopted a simple cut-based method in 
this analysis. One would instead invoke more complex discriminators to enhance the significances within the designed 
luminosity, for example, multi-variate analyses. 

To conclude, after the first few years of the LHeC running, by adopting more complex discriminator and/or multi¬ 
variate analyses, we expect that both h and H signals will appear at the LHeC. 
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